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ABSTRACT 

Distributions and linear regressions are discussed. 
The section dealing with the former topic emphasizes the usefulness 
of computer programs in statistics, demonstrating their ability to 
handle tedious and time--consuming tasks* The normal curve is stressed 
since the assumption of a normal distribution is common. An example 
of 200 data points is presented which illustrates a computer 
program's ability to do descriptive statistics for data not grouped, 
data which is distributed into any number of classes of equal width 
and testing if the 200 points are normally distributed. The section 
of linear regression emphasizes how computer programs can lead to 
erroneous results. In the area linear regression there are in current 
use "canned programs" which display the t-values corresponding to the 
coefficients of the least square estimate, but do not take 
autocorrelation into account. In the presence of autocorrelation it 
is incorrect to use these t-values. An example is presented in which 
autocorrelation is present but can be removed by a suitable 
transformation of variables, (Author) 
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This paper is divided into two sections, DISTRIBUTIONS, and LINEAR REGRESSION. The 
first section DISTRIBUTIONS emphasizes the usefulness of computer programs in statistics. 
Some very tedious and time-consuoiiiig tasks are made considerably easier and clarified by use 
of these computer programs. In particular, the normal curve is emphasized since in many 
areas of statistics the assumption "a population is normally distributed" appears time and 
time again. An example of 200 data points is presented illustrating the computer program's 
ability to do descriptive statistics for data not grouped, data which is grouped into any 
number of classes of equal width, and testing if the 200 points are normally distributed. A 
graph of the data also appears in the output. 

The second section LINEAR REGRESSION emphasizes tow computer prograuis can lead to 
erroneous results. In the area of linear regression there are in current use ''canned 
programs" which display the t-values corresponding to the coefficients of the least square 
estimate, but do not take autocorrelation into account. In the presence of autocorrelation 
it is incorrect to use these t-values. An example of this is presented where autocorrelation 
is present but can be removed by a suitable transformation of variables. 



Distributions 

The following are the salient features of a program written to obtain the following: 

1. Descriptive statistics for data which is not-grouped. 

2. Indicators that a set of points are normally distributed (necessary but not 
sufficient conditions) . 

3. Grouping original data into any number of classes of equal width, with descriptive 
statistics. 

A graph of the grouped data (frequency plotted against mid-point of CLASS) . 
5. CHI-SQUARE test for goodness of fit of normal curve. 
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The following is a display of an input and output cf the program for data points as 
indicated : 
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CHOOSE WHAT INFORMATICN YOU NEED 

A) NOT-GKOUPED EAXA 

B) GROUPED DATA 

C) BOTH 

D) NONE OF THE ABOVE 

ANSWER A,B,C, OR D 
?C 

HOW MANY CLASSES THERE? 
?9 

WHAT IS THE MINIMUM NUMBER OF THE FIRST CLASS? 
?3.67 

WHAT IS THE MAXIHUM NUMBER OF THE LAST CLASS? 
?4-75 



STATISTICS FOR •NOT-GROUPED' EATA 



THERE WERE 200 OBSERVATIONS MADE 

THE AVERAGE IS 4.1714 

THE VARIANCE IS 3.09021E-02 

THE STANDARD DEVIATION IS .17579 

THE MAXIMUM VALUE IS 4.7 

THE MINIMUM VALUE IS 3.68 

THE RANGE OF THE OBSERVATIONS IS 1.02 

THE MEASURE OF SYMMETRY IS .260124 

THE MEASURE OF PEAKEDKESS IS 3.10359 



THE FIRST DEVIATION BAND BETWEEN 3-99561 AND 4.34719 
CONTAINS 68% OF THE OBSERVATIONS 

THE SECOND DEVIATION BAND BETWEEN 3.81962 AND 4.52298 
CONTAINS 95% OF THE OBSERVATIONS 

THE THIRD DEVIATION BAND BETWEEN 3*64403 AND 4.69877 
CONTAINS 99. 5« OF THE OBSERVATIONS 
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STATISTICS FOR 'GEJUPED* DATA 



CLASS BOUNDARIES 

FROM TO MIDPOINT FREQUENCY 

3.67 3.79 3.73 3 

3.79 3.91 3.85 9 

3.91 4.03 3.97 28 

U. 03 a. 15 4.09 54 

4. 15 4. 27 4.21 51 

4.27 4.39 4.33 31 

4.39 4.51 4.45 17 

4.51 4.63 4.57 6 

4.63 4.75 4.69 1 
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3.73 + * 

3.85 ♦ * 

3.97 ♦ * 

4.09 ♦ * 

4.21 ♦ * 

4.33 ♦ * 

4.45 ♦ * 

4.57 + 

4.69 +♦ 



+++++++++++++++++++++++++++++++++++++++++++++++++++++ 

THE NUMBER OF OBSERVATIONS MADE WERE 200 

THE AVERAGE IS 4. 1728 

THE MODE IS 4.09 

THE VARIANCE IS 3.20242E-02 

THE STANDARD DEVIATION IS .178953 

THE MEDIAN IS 4. 16412 
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36. 108 
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17 


15.958 


6. 80388E-02 
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4.498 


. 501557 


1 


.8088 


4. 51996E-02 



TOTAL 



2.41273 (Ccmputed Chi-square value) 



1. The usual measures such as average, variance, and standard deviation are obtained- 
In addition the maximum value, minimum value and range are determined. 

2. Preliminary indications are that the set of 2C0 numbers is normally distributed 
since the measure of sytDmetry is approximately zero, peakedness is approximately 3, 
and one, two, and three standard deviation bands contain respectively 68%, 95%, and 
99. 5% of the data. 

3. Nine classes were asked for and the grouped data for these nine classes is 
presented with the midpoint of each class. Ihe usual measures such as average, 
variance, standard deviation, median and mode are given for the data grouped into 
these nine classes. It is noted that the average and standard deviation for the 
grouped data {U.1728 and .178953) are approximately the same for the data when not- 
grouped {a.171£4and . 17579). 

4. The graph indicates the data is approximately symmetric about the average value of 

a. 1 71 a. 

5. The CHI-SQDARE value of 2.41 273 indicates th^t the normal curve (average of U. 11728 
and standard deviation cf .176953) fit is good at both the 5 percent and 1 percent 
level of significance { 0(1) is observed frequency and E(I) is expected frequency)- 



The program is particularly valuable in that it will offer descriptive statistics for 
any number of classes asked for. In seme realistic situations it is important to compare 
results using different numbers of classes and this can new be done with considerable ease. 
The graph of the data is of particular importance in that ore can quickly get an initial 
impression of how the data behaves. The program emphasizes the normal distribution. The 
reason for this is that in such topics as linear regression (second section of this paper) an 
assuaiption often stated is that a set of points is normally distributed- 



An assumption in the linear regression model is that the error terms are statistically 
independent. If violated, it is referred to as a problem cf autocorrelation, invalidating 
use of the usual t-tests. There are available computer "canned programs" currently which 
display the t-values corresponding to the coefficients of the least square estimate, but do 
not take autocorrelation into account. Therefore, in the presence of autocorrelation it is 
incorrect to use these t-values [2, p. 80]. 

Under suitable conditions, depending on the given data, autocorrelation can be 
eliminated by transformation to a new set of variables. In the transformed system it is 
valid to use the t-values. This subject is treated clearly by Frank [3]. A typical case 
where autocorrelation is not considered is Parker and Segura's article [^], an interesting 
and useful paper. In their article, sales are forecast in the home furnishing industry using 
new mar ri ages durin g the year , housing starts, annual disposable income, and time trend as 
explanatory variables, using 2^4 years of data. This is an example of time-series data where 
one must be especially careful for the presence of autocorrelation. In the transformed 
system it is demonstrated that autocorrelation is eliminated and the least square estimate 
provides a reasonably good estimate of the original set of data with the t-values significant 
at approximately a 5% level of sigr.if icance. 
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GIVEN 2^ YEARS OF DATA 



Disposable 
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1965 
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Define the fcllowirg for year i: 

Hj^ = Housing Starts 
I^=Annual Disposable Income 
M . =New Marriages 
S^=Gross SaleE 
Ti =Tini€ Trend 
i*-1,2,...,24 

Presence of autocorrelaticn in error_ternts 

Using a linear regression computer program the fcllcwing least square estimate is 
obtained: 

S = 50.605 + .036H + 1.2211 -.068M -19.483T (1) 
The following table is then established: 



GROSS SALES 



ESTIflATED 
GEOSS'SALES 



ERROR TERMS 
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92. 94 0 


95.673 


-2.752 
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122.440 


116.729 


5.711 
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125. 57C 


135.752 


-10. 181 
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11 0. 460 


134.467 


-24.006 
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139.400 


146.578 


-7. 178 


6 


154. 02C 


143. 129 


10.891 
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133.008 


6.122 


10 


156.330 
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16.824 


11 


140. 470 


130.447 


10.024 


12 


128.240 


131. 561 


-3.341 


13 


117.450 


135.877 


-18.427 
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143. 38^J 


-10.744 


15 
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127. 395 


-1.236 


16 


116.990 


126.037 
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17 


123. 900 


134.937 


-11 .037 


18 


141.320 


139.232 


2 .088 


D9 


156.710 


151 .599 


5.111 


20 


171 .930 


165.210 


6.720 


21 


184.790 


174.802 


9.988 


22 


202.700 


205. 204 


-2 .504 


23 


237. 340 


236.319 


-.979 


2a 


254.930 


259. 474 


-4. 544 



The t-values corresponding to the coefficients of the least square estimate are not displayed 
here since it will be shown that there is the presence of first order autocorrelation 
invalidating the use of these t-values. The main point cf this section is that a computer 
program is incorrect if it displays these t-values if autoccr relaticn is present. 

Denoting by Xi the error term for year i , a test for first order autocorrelation is 
based on the Durbin-Hat son statistic [3, p. 276]. 




Since the number cf observations is 24 and 4 independent variables are under 
consideration it follows frcm Table E of Frank's book [3, p. 276] that the d^ and ^2 for a 
two-tailed critical region of 5 per cent are given by 



d^ = .91 d2 1.66 (3) 



and that therefore 



d < dj^ < d2 (4) 
roncluding that there is significant first-order autoccr relati on [3, p. 279]. 



Transformation to. d new__set of variables--eliiEinating autoccr relation 

The first step is to determine ^, the coefficient cf autcccrrelation. The estimated 
value of p is given by Frank [3, pp. 280- 281]. 

24_ ^ 
^ ^2 ^ ^'^ 

P ^ ~^ ^ (^J 



An approximation to ^ is given by 

^ = 1/2 (2-d) = .557U (6) 

where d is the Durbin- Wat son statistic. 

Denoting by Hj^, 1^, H^, S-^ , a new set of variables the transformation is given by 
[Frank [3, p. 280} 



H. = 


«i 


-. 5571) 




= 


li 


-. 5571) 


Ii-1 


«i = 


"i 


-. 5571) 


«i-l 


Si = 


Si 


-. 557H 


^i-l 


T. = 


■^i 


-. 5571) 


^i-l 



where i = 2,3, . . . ,24. 



Performing these transformations 


leads to the 


following set 


of columns: 
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H. 
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\ 
I . 




M . 
1 


s. 

1 




T? 

1 


1 

2 


527.303 


80.93 1 




7 14. 022 


7C.647 




1 . 443 


3 


507.940 


93. 823 




701. 239 


57.323 




1.885 


a 


562.217 


82.244 




570. 568 


40.46S 




2. 328 


5 


914.691 


101. 457 




786. 325 


7 7, 831 




2.770 


6 


347.604 


110. 189 




665. 833 


76.320 




3.213 


7 


576. 002 


109. 643 




649. 964 


71.741 




3. 656 


8 


527.706 


116. 91 4 




688. 178 


64.391 




4.098 


9 


671 . 198 


116.223 




628. 277 


54. 279 




4. 541 


10 


720.295 


132. 600 




700. 490 


7 8.781 




4. 984 


11 


417.851 


139. 952 




731. 637 


53.333 




5* 426 


12 


465.836 


145. 241 




634. 538 


4S. 944 




5. 869 


13 


696.101 


146. 546 




604. 883 


45. 970 




6.311 


14 


783.541 


160. 172 




685. 229 


67. 175 




6.754 


15 


420. 637 


161.770 




694. 261 


52.228 




7. 197 


16 


616.886 


169.314 




695. 867 


46.670 




7,639 


1" 


729.345 


182. 187 




717.719 


58.691 




8,082 


18 


796. 196 


189.838 




773. 325 


72.260 




8. 524 


19 


637.817 


211 .080 




797. 079 


77.94C 




8.967 


20 


630.736 


225. 744 




830. 848 


84.582 




9. 410 


21 


343.605 


243. 829 




846. 831 


8e.958 




9.852 


22 


645. 183 


264. 65 1 




835. 175 


9S.700 




10.295 


23 


799.952 


286. 398 




992. 715 


124.357 




10.737 


2a 


629.539 


300. 23 9 




984. 336 


122.639 




11. 180 


For 


these transformed 


varia blefs 


the 


least square estimate is 


given 


by 




S = -43. 879 


+ .023H + 


. 7131^ + . 088H^ 


-12.709T^ 


(8) 




In 


the transformed system the error 


terms are 


given as fellows : 
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i 



(ESTIMATED) 



ERROR TERMS 



1 




— 




2 


70.647 


70.590 


.0 57 


3 


57.323 


72. 583 


-15,260 


^ 


40.469 


48.479 


-8.010 


5 


77 .83'' 


83.754 


-5.923 


6 


76.320 


60.525 


15.795 


7 


71 .741 


58. 447 


13,294 


8 


64 .391 


60. 240 


4. 152 


9 


54.279 


52. 205 


2. 074 


10 


78.781 


65.755 


13.0 26 


1 1 


53.333 


61.051 


-7,713 


1 2 


49 .944 


51.782 


-1 .838 


13 


45 .970 


49.856 


-3.885 


ia 


67.175 


63.053 


4 . 121 


15 


52.228 


50. 890 


1 .337 


16 


46.670 


55.368 


-8.698 


17 


58.691 


63. 470 


-4 . 779 


18 


72.260 


69.751 


2 . 508 


19 


77.940 


77.667 


. 273 


20 


84 .582 


85.304 


- .722 


21 


88.958 


87. 279 


1 .679 


22 


99.700 


106.915 


-7,215 


23 


J 42 .357 


1 29.868 


12 . 489 


2^ 


122.639 


129.398 


-6 .759 



For these error terms the Durbin-Watson statistic is given by 

d = 1.8015 (9) 

Since 23 observations and 4 independent variables are under consideration one 
establishes from Table E cf Frank's book [3, p. 276 ] that 

d^ = .89 a 2 = 1-67 (10) 

And therefore in the transformed system 

d > d 2>dj^ (11) 

in which case one accepts the hypothesis cf no first crder autocorrelation [3, p. 279]. 
Therefore in the transformed system it is valid to use the t-values. 



Siflnif ica nee _of_t- values 

The transformed equation with the corresponding t-valu€s below them are 

S = -43.879 f .023H ♦ 7l3l ♦ 088M -12.709T 

- 1.952 1. 694 2. 959. 1 . 658 3. 237 (1 2) 

Coefficient cf Determination = , 893 

Standard Error of Estimate = 8.845 

In order to determine the significance of the t-valuos one uses a Student's-t 
disbributicn with (N-K) degrees of freedom. Since in the transformed system 

N = number of observations = 23 

K = total number of variables under consideration = 5 

It follows that the number of degrees of freedom is 18- The critical value for a 5% 
one-tailed test of significance using the Student' s-t distribution with 18 degrees of freedom 
is 1.729. The t-values associated with the constant term, dispcsable personal income and 
time trend are significant at a 5^ level. The t-values associated with new housing starts 
and new marriages are significant at approximately a S% level. 



NOTES AND REFERENCES 

1. Ya-Lun Chou, Statistical Analysis. Holt, Rinehart, and Winston, Inc., 1969. 

2. w. L. Hays and L. W. Winkler, Statistics, Vol. II. Holt, Rinehart and Winston, Inc., 
1970. 

3. C. R. F. Frank, Jr., Statistics and Econometrics. Holt, Rinehart and Winston, Inc., 
1971. 

M. G. C. Parker and E. L. Segura, How to get a tetter £or*5cast. Harvard Business Review, 
1971, March-April, 99-109. 

♦Please send correspondence to Dr. Munick. 



ERIC 



